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^ ■ Abstract 

We calculate electric and magnetic masses of gluons between T — Tc 
\^ ' and QTc using lattice QCD (quantum chromodynamics) in the quench ap- 

, proximation for color 5f/(3). We find that magnetic mass has finite val- 

ues in this region, and the temperature dependence of the electric mass 
' is consistent with that determined using the hard-thermal-loop pertur- 

^ , bation. The hard-thermal-loop resummation improves significantly the 

magnitude of the electric mass comparing to the leading order perturba- 
I tion. Both screening masses have little gauge parameter dependence. 

OO , We are expecting relativistic heavy-ion coUisions to soon yield evidence of a 

■ quark-gluon plasma (QGP), by creating sufficiently high temperature to bring 

the system into the deconfining phase. This new form of matter is a prediction 
of QCD and will be studied for the first time in the laboratory. One of the 
\ most urgent theoretical tasks is to investigate the characteristics of gluons, the 

I ' component of QGP that controls the force, at a finite temperature based on 

Qh! QCD. If we can calculate these features in a reliable manner, phenomenolog- 

ical observations, such as the jet quenching, may be analysed based on these 
knowledges. 

^ ' There has been essential progress in the perturbative method of treating 

Jv>( \ hot QCD. Using the hard-thermal-loop resummation technique[0, we can now 

^ ■ handle plasma effects of gauge theories. The correct magnitude of the free 

energy from lattice QCD was reproduced by such a calculation^, H|. 

Fundamental quantities of the hot gluon plasma have been studied pertur- 
batively. Electric mass is given by the lowest perturbative calculation as 

mefl = \I^9T- (1) 

Magnetic mass, m^, is, however, beyond the scope of the perturbative treatment| 
Ij. It vanishes in the lowest perturbation, but if indeed to„ = 0, the gluons 
are essentially unscreened, even if the electric mass is finite Also, zero mag- 
netic mass causes serious infrared divergence and makes consistent theoretical 
treatment very difficult (see Chap. 10.2 of Ref.Q). 

The maximum temperature accessible at RHIC or LHC is of the order of 
IGeV or 5Tc which is far below the regions where the perturbative calculation 
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can be safely applied. Therefore nonperturbative analysis is highly desirable in 
order to gain a quantitative understanding of the properties of QGP. Since the 
pioneering work for temperature dependencies of screening masses by Heller, 
Karsch and Rank for SU(2) color QCD|8|, there have been several activities, 
but no SU{3) lattice simulations of finite temperature gluon propagators in the 
literature. Indeed, calculations of gauge-dependent observables are very time- 
consuming and often unstable. 

In this letter, we report the electric and magnetic screening masses of gluons 
at finite temperature measured in lattice simulations based on the SU(3) gauge 
theory, i.e., QCD in the quench approximation. We calculate gluon propagators; 
which allows direct comparison with perturbative calculations. 

We calculate gluon correlation functions with finite momentum, 

G^,{PuP^,Py,z) = (TrA^(P„Py,P(,z)A,(-P„-Py,-Pt,0)), (2) 

where gauge potentials ^^'s are evaluated from the link variables U^{x) G SU(3) 
as A^{x) = {Ufi{x) — Uj^{x))/2gi. The calculation of propagators with finite mo- 
menta enables us to construct transverse and longitudinal components without 
ambiguity Q; we define electric and magnetic parts of the gluon propagators as 
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G„(|^,0,0,z)-f Gtt (0,^,0, z) 

Gxxi^l "^J Oj ^) + C!yy{jf-, 0, 0, Z) 
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which correspond to standard projected propagators in perturbative calcula- 
tions. Here P^'s are lattice momenta, and the size of the lattice is given by 
NxNyNzNt- At a sufficiently large distance, they are expected to behave as 

Ge(„)(^)cxe-^=<'")^ (4) 

In order to quantize the theory and to extract gauge-dependent observables, 
we employ stochastic quantization a la ZwanzigerJlOj, 

dAf, 5S 1 , K 

where r is the Langevin time, the second term of the r.h.s. is the gauge fixing 
term, and rj represents Gaussian random noise. The gauge parameter is given 
by a, and a = corresponds to the Lorentz gauge. The lattice version of the 
algorithm was proposed in Ref . as 

f7^(a;,r-|- Ar) ^ uj'' {x,T)e^^^^^'^ Uf,{x,T)uj{x + ji.r), (6) 

Here, lo means a gauge rotation matrix and / a driving force. 

The gauge rotation step and Langevin step are executed alternately. 
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There are two reasons for using the stochastic gauge fixing method here: 
one is practical and the other is conceptual. Contrary to the standard Wilson- 
Mandula gauge fixing method where an iterative procedure is applied and 
the number of iterations is unpredictable for large lattices, in algorithm (|^), we 
repeat the steps of update and gauge rotation one after the other. There is 
no convergence problem of gauge fixing. Therefore we can precisely estimate 
CPU time. Moreover, this algorithm, which may change the gauge parameter a, 
has an advantage in testing gauge independence. A conceptual problem is the 
Gribov ambiguity appearing in non-Abelian gauge theories Although our 
algorithm may not eliminate copies completely, Zwanziger's stochastic gauge 
fixing algorithm maintains the gauge system inside the Gribov region[|l0[ ^ if 
we start from the trivial configuration {A^ — 0}. 

We perform simulations on an N^NyNzNt = 20 x 20 x 32 x 6 lattice with 
the standard Wilson gauge action. The physical temperature is tuned by chang- 
ing the lattice cut-off scale l/a[GeV], i.e., gauge coupling. The relation between 
gauge coupling and the lattice scale is estimated using the data in Ref . , which 
covers the range from /? = 6/5^ = 6 to 7.1. The confinement/deconfinement 
transition temperature is set to be 256 MeV[p^. Then we study the regions 
of T/Tc = 1 ~ 6. The Runge-Kutta algorithm is applied to reduce the finite 
Langevin step At dependence [Q. We perform the simulation for a set of pa- 
rameters with At = 0.03 ~ 0.05 and extrapolate all results linearly to At ~ 0. 
The gauge parameter a is fixed to be one, except in the gauge dependence test. 

We observe a large fluctuation of gauge propagators, particularly at large 
distances. We show in Fig|| a typical behavior of gluon propagators G{z) as 
a function of the Langevin step. In order to estimate the necessary number of 
Langevin update steps, we wait until < Gxx{z) > and < Gyy{z) > calculated in 
Eq.H becomes equal to each other. A typical number of simulations is 0.2 0.4 
million steps. 

The dynamics of a gluon propagator is completely different in confinement 
and deconfinement regions In Fig.^j, the gluon behavior is massive and we 
see that as distance z increases, the mass of the propagator in confinement 
regions seems to become infinite; the gluon is completely screened in the con- 
finement phase. Consequently, we cannot employ the assumption, Eq.(Q). On 
the contrary, the propagator decreases exponentially at a long distance with a 
finite damping rate in the deconfinement region. 

To extract screening masses from the propagators we use the following fitting 
functions under the lattice periodic condition: 

Ge(™) ~ cosh(Se(„)(p)(z - NJ1)). (8) 

We employ values for z > 1/T{= Nta), because a screening effect occurs in a 
sufficiently long range. Then, we obtain static mass from E(jj). 

To get the final result at At = 0, we must extrapolate with respect to 
Langevin step width. Fig.^ represents E{p) measured here versus At. The 
slight dependences of At cause us to use a linear function when fitting data. 

Screening masses are physical and expected to be gauge independent. How- 
ever, since the definition of the gluon propagator, Eq.(^ itself, is generally gauge 
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Figure 1: Typical gluon behavior with Langevin steps. We find that fluctua- 
tions of Gxx{^'2) are much larger than those of Gxx{Q)- In order to investigate 
the screening effect, long-range contribution should be adopted. Consequently 
we need a large number of statistics. 



dependent, it is nontrivial whether the screening masses are gauge invariant or 
not. In addition, since we fail to define the magnetic mass by a perturbative 
analysis, it is particularly important to check its gauge invariance. In Fig.^, we 
show the gauge parameter a dependence of electric and magnetic masses. The 
result strongly suggests that they are gauge independent and physical observ- 
ables. 

We study the temperature dependence of screening masses at T/Tc = 1 ~ 
6 which would be expected to be realized in high-energy heavy ion collision 
experiments such as RHIC or LHC. Fig. 3 shows the behavior of electric and 
magnetic masses as a function of temperature. The magnetic part definitely 
has nonzero mass in this temperature region. As T increases, both masses 
decrease monotonically, and at almost all temperatures, magnetic mass is less 
than electric mass, except very near Tc where the electric mass decreases very 
quickly as T approaches Tc- 

We fit the data above T - LST^ by 

^ (9) 
which are predicted through the perturbative and 3-D reduction analysis H- 
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Figure 2: The electric propagator behavior in confinement (circle) and decon- 
finement phase (the others) with = j^. We find all propagators become 
massive comparing with the free propagator (long dashed line). In confine- 
ment regions, the propagator at long distance seems to have infinite mass and 
vanishes, whereas we find that the propagator beyond Tc, has finite mass. 



Here, we use running couplings as 
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bi log(21ogA^/A) \ 
2bo log(Ai/A) J ' 
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and we set /i — 2ttT which is a Matsubara frequency, as the renormaliza- 
tion point and A = l.OSTcQ as the QCD mass scale, bo — llNc/^Sir'^ and 
bi = (34/3)(A*'c/(167r^))^ are the first two universal coefficients of the renormal- 
ization group. We obtain Ce = 1.63(3), xVNDF = 0.715 and C„ = 0.482(31), 
X^/NDF = 0.979. The scaling expected in Eq.(^ works well for the electric 
mass although the magnitude Ce is larger than the leading order perturbation 
Eq.(0). 

The hard-thcrmal-loop rcsummation technique applying the free energy of 
hot gluon plasma has been successful, i.e., it gives the correct sign and roughly 
the correct mag nitude for T > 2rc[|. Rebhan gave a formula for the electric 
mass in the one-loop resummed perturbation theory p9|. 
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Here, we assume the magnetic mass to be of the order of g^. Substituting our 
fitted value for m„j, we can solve the above equation iteratively. In Fig. 3, we 
show this hard-thermal-loop resummation result together with the lowest order 
calculation. The hard-thcrmal-loop result gives a better description than naive 
perturbation, upon comparing our numerical experiment. 
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Figure 3: At dependence of masses is slight. Therefore, obtaining a final value 
at At — 0, we use the linear function for extrapolation. Here, it is important 
to apply the Runge-Kutta algorithm in order to suppress At error. 



In Ref . pO| , the magnetic mass was estimated as rUg = 2.38A^c5^r/47r = 
0.568(7^T using a self-consistent inclusion technique. This is very close to our 
fitted result Cm = 0.515(35). 

The electric mass was obtained also from Polyakov line correlation functions 
at finite temperature |22j and those contradict our results. 3-D reduction 
argument|^^ has shown that rrie/gT goes down when T increases, but even at 
T ~ 1000A-|y-gr the electric mass is still about 3TOe,o- 

In conclusion, we have measured gluon propagators and obtained the electric 
and magnetic masses by lattice QCD simulations in the quench approximation 
for SU(3) between T = Tc and 6Tc. Features of QGP in this temperature region 
will be extensively studied theoretically and experimentally in near future. We 
observed that the magnetic mass does not vanish there; To our knowledge, this 
is the first reliable measurement in SU(3). It can be approximately fitted by 
Cmg^T, and that the electric mass is consistent with the hard-thermal-loop 
resummation calculation. We also confirmed that both electric and magnetic 
masses are gauge independent. The calculational technique presented here can 
be applied to the study of quark screening effects at finite temperature. 

We would like to thank A. Niegawa, S. Muroya and T. Inagaki for many 
helpful discussions. Numerical calculations were carried out on SX5 at RCNP, 
Osaka Univ., VPP5000 at SIPC, Tsukuba Univ. and HPC computer at IN- 
SAM, Hiroshima Univ.. This work is supported by Grants- in- Aid for Scientific 
Research from Ministry of Education, Culture, Sports, Science and Technology, 
Japan (No.ll694085, No.ll740159, and No. 12554008). 
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Figure 4: Gauge dependence graph for electric and magnetic screening masses. 
Gauge dependences of both screening masses are very shght in nonperturbative 
regions. 




Figure 5: Temperature dependences of electric and magnetic screening masses. 
Dotted line was fitted based on the scaling of nig ~ g'^T. For the electric 
mass, the dashed and solid lines represent a leading order perturbation and the 
hard-thermal- loop resummation result, respectively. 
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